'  AD-A126  222 
UNCLASSIFIED 


fi  COMPARISON  OF  OPTICAL  VERSUS  HARDWARE  FOURIER  1/ 

TRANSFORMS(U)  VIRGINIA  POLVTECHNIC  INST  AND  STATE  UNIV 
BLACKSBURG  DEPT  OF  PHVSICS  S  P  ALMEIDA  21  OCT  82 
AFOSR-TR-83-1077  AFOSR-83-0200  F/G  12/1  NL 


Submitted  to 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
Final  Technical  Report 
(Grant  No.  AF0SR-83-0200) 

October  31,  1983 

A  COMPARISON  OF  OPTICAL  VERSUS  HARDWARE 
FOURIER  TRANSFORMS 

by 

Silverio  P.  Almeida 


UNCLASSIFIED 


SECURITY  classification  of  This  PAGE  (Whan  Data  Enteted) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3  RECIPIENT’S  CATALOG  NUMBER 


♦  TITLE  (and  Subtttla) 


A  comparison  of  optical  versus  hardware  Fourier 
trans  forms 


B 


A'j’’-'OR(»i 

*.  Silverio  P.  Almeida 


6.  PERFORMING  o^g.  report  number 


a.  contract  or  grant  numbers; 


AFOSR-83-0200 


}  pra^ORMING  ORGANIZATION  NAME  ANO  ADDRESS  10  PROGRAM  ELEMENT.  PROJECT,  TASK 

...  .  „  „  .  ,  _  .  ,  _  ...  .  AREA*  WORjt  UNIT  NUMBERS 

.  lr  gima  Polytechnic  Institute  and  State  University  L/IOr^r 
Contract  and  Grant  Administration,  301  Burruss  Hail  2305/D9 
Blacksburg,  VA  24061 


I!  CONTROLLING  OFFICE  NAME  AND  ADORESS  '2-  REPORT  DATE 

.-.tOSR/ARB/^es  October  31,  1983 

Building  410 

Bolling  AFB,  DC  20322  _ 


1  4  MONITORING  AGENCY  NAME  *  ADDRESS!-//  different  from  Controlling  Otllca)  IS.  SECURITY  CLASS,  (ol  thie  report; 


AFOR/NE 

Building  410 

Bolling  AFB,  DC  20332 


Unclassified 


IS*.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16  Distribution  STATEMENT  (ol  Ihla  Report) 


Approved  for  public  release? 
distribution  unlimited* 


IT  DISTRIBUTION  STATEMENT  (ol  the  abatracl  entered  In  Block  20,  II  dlllerenl  from  Report) 


'8  SUPPLEMENTARY  NOTES 


AFOSR  Final  Report 


19  <Ev  #ORDS  (Continue  on  ravaraa  aide  if  nacaaamry  and  idantlfy  by  block  numbar) 

Optical  Fourier  Transforms 

Array  Processing  of  Fast  Fourier  Transform 

Spatial  Filtering 


20  ABSTRACT  f Contlnua  on  ravaraa  aide  If  nacaaaary  and  idantlfy  bv  block  numbar) 


A  comparison  of  the  optical  transform  versus  the  hardware  Fourier  transform 
was  made  for  two  input  functions:  a  rectangular  aperture  and  a  six-pointed 
star  design.  The  optical  transforms  were  modified  using  sharp  and  smooth  lev- 
pass  as  well  as  highpass  filters.  The  results  were  compared  against  those 
obtained  using  the  array  processor  hardware  at  the  Brooks  Air  Force  Base, 
School  of  Aerospace  Medicine.  The  "ringing"  observed  in  the  lowpass  filtered 
results  of  the  array  processor  was  simulated  by  the  optical  methods,  However, 
certain  asymmetries  in  the  hardware  "ringing"  data  were  not  observed  optically. 


'  1473  EDi*lON  3F  I  MOV  83  IS  OBSOLETE 


ABSTRACT 


A  comparison  of  the  optical  transform  versus  the  hardware  Fourier 
transform  was  made  for  two  input  functions:  a  rectangular  aperture  and 
a  six-pointed  star  design.  The  optical  transforms  were  modified  using 
sharp  and  smooth  lowpass  as  well  as  highpass  filters.  The  results  were 
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the  Brooks  Air  Force  Base,  School  of  Aerospace  Medicine.  The  '"ringing” 


observed  in  the  lowpass  filtered  results  of  the  array  processor  was 


simulated  by  the  optical  methods.  However,  certain  asymmetries  in  the 
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hardware  "ringing*  data  were  not  observed  optically.  A  discussion  of  the 


ringing”  and  how  to  eliminate  it  is  also  presented. 
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I.  INTRODUCTION 


The  main  purpose  of  this  research  is  to  compare  the  Fourier  transform 
and  its  inverse  filtered  Fourier  transform  obtained  with  the  Digital  Image 
Processing  (DIP)  hardware  system  located  at  the  School  of  Aerospace 
Medicine,  Brooks  Air  Force  Base,  San  Antonio,  Texas  against  the  optical 
version  of  the  transforms  obtained  at  the  Virginia  Tech  coherent  optics 
laboratory.  The  Air  Force  test  patterns  consisted  of  a  rectangular  aper¬ 
ture  and  a  six-pointed  star  design.  They  were  both  recorded  on  35  mm 
transparencies,  and  provided  to  us  by  Dr.  Ralph  G.  Allen,  Director  of 
the  Laser  Effects  Branch  (Division  of  Radiation  Sciences).  The  DIP 
system  consisted  of:  an  Analogic  array  processor  on-line  to  an  Aydin 
monitor  and  a  PDP-11-32  mini-computer  system.  Both  the  rectangular 
aperture  and  six-pointed  star  design  were  digitized  (512x512)  pixels  and 
entered  into  the  Analogic  array  processor  which  performed  the  hardware 
Fourier  transform  (HFT).  This  HFT  was  then  operated  on  with  lowpass  and 
highpass  Modulation  Transfer  Functions  (MTF) ;  next  the  inverse  hardware 
Fourier  transform  (IHFT)  was  performed.  The  Aydin  monitoT  was  used  to 
display  the  HFT  and  the  IHFT.  The  results  were  photographs  of  the  Aydin 
monitor  screen  displays  and  recorded  on  35  mm  transparencies.  It  is  this 
set  of  transparencies  including:  the  original  input  .designs  (rectanuglar 
aperture  and  six-pointed  star);  the  HFT  and  IHFT  that  we  wished  to  compare 
against  their  optical  equivalents.  In  some  cases  we  were  able  to  duplicate 
the  hardware  results  using  optical  methods.  In  those  instances  where 
differences  in  the  two  methods  were  observed  we  have  recorded  the 
differences  with  explanations  for  them. 
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It  is  important  to  test  the  DIP  system  Fourier  transform  capabilities 
against  the  optical  (traditional)  methods.  The  DIP  system  could  be  very 
useful  for  on-line  filtering  (Modulation  Transfer  Manipulation)  of  Images 
and  rapid  analysis. 


II.  OPTICAL  FOURIER  TRANSFORMS :  OFT 

A  two  dimensional  Fourier  transorm  can  be  performed  using  a  convergent 

lens  and  coherent  illumination  of  the  input  pattern^.  Various  optical 

set-ups  can  be  implemented  to  carry  our  the  OFT.  For  this  experiment  we 

2 

have  chosen  a  converging  beam  approach  rather  than  the  classical  2-f 
3 

configuration  .  Both  methods  are  equivalent  for  our  purposes  as  is  shown 
in  ref.  (2).  Figures  (1,  2)  show  the  converging  beam  schematic  used  to 
obtain  the  OFT.  The  high  and  lowpass  spatial  frequency  filters  were 
located  on  glass  microscope  slides  as  either  a  black  dot  (highpass)  or  a 
black  annulus  (lowpass)  and  positioned  in  the  Fourier  transform  plane  via 
micrometers  with  X-Y  movements.  Since  the  exact  MTF  cut-off  frequencies 
used  in  the  IHFT  are  unknown  for  the  transparencies  provided  to  us,  we 
have  tried  various  optical  pass  filters  in  order  to  simulate  these  results. 


A.  Rectangular  aperture  Fourier  transform  theory 

1  2 

The  analytic  treatment  of  the  OFT  can  be  outlined  as  follows  '  . 
Given  a  rectangular  aperture  whose  amplitude  transmittance  is  given  by 

x  y 

t(x  yQ)  -  rect  (j2-)  rect  (^-) 
x  y 


(1) 


2- 


where 


rect  (x  ) 
o 

and  similarly  for  rect  (yQ) •  i*y  represent  the  aperture  widths  in  the 

x  and  y  directions,  respectively.  Assume  that  the  pupil  function 
o  o 

2  2  2 

p(x  ,  y  )  ■  1  since  x  +  y  <  a  where  a  *  radius  of  the  lens;  also  that 
o  Jo  o  o 

the  amplitude  illumination  is  unity  (i.  e.  A«l).  Then  in  the  Fresnel 
approximation  the  field  distribution  in  the  Fourier  transform  plane 
g2(xf,  yf)  is  given  by: 

A  exp  [i  +  fy>]  ,fs 

*2  (lf’  y£>  *  — : — ns -  '5'  • 


where  A  =  1  (amplitude  transmittance), 

PHI  (pupil  function) 

f  m  d  (input  plane  is  assumed  to  be  next  to  the  backside  of  the  lens; 

f  *  d).  f  is  the  focal  length  of  the  lens. 

The  Fourier  transform  of  the  transmittance  function  t(xo,yo)  is  given 
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F[t(v  yo)l 


L  L  sine  (L  f  )  sine  (L  f  ) 

Xy  XX  y  y 


“h“e:  E  <xf>  •  fy  5  <!?> 


are  Che  spatial  frequencies,  some  of  which  will  be  passed  via  the 
previously  mentioned  low  and  highpass  filters.  We  now  get  that  the 
Fourier  transform  of  the  rectangular  aperture  is  given  by 


g2(xf*  yf5 


exp[i  ^  <x*  +  y*)] 


L  L  sine  (l  f  )  sine  (L  f  )  (4) 

xy  x  X  yy 


The  photograph  representing  this  Fourier  transform  is,  however,  a 
square-law  detector.  That  is, it  records  the  intensity  distribution  of 
(4)  given  by: 


l?  l}  ,  Lx  2  L  y 

■Sleet  (xf*  yf5  *  x2f2  sinc  (Tf_)  3inC  (  Xf  J 


B.  Star  hperture  Fourier  transform  theory 


The  six-pointed  star  can  be  treated  as  a  two  dimensional  triangle 
function.  That  is,  by  using  the  following  Fourier  transform  properties: 


(a)  Rotation:  f(r,  9  +  0q)  -*-*•  F(ui,  $  +  0^) 


(b)  Distributive:  F  {  f^x.y)  +  f2(x,y)}»  FCf^x.y)}  +  F{f2(x,y)} 
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(c)  Scaling:  f(ax,  by)  1  v  , _ x_  , 

[ST  (  a 


(d)  Translation:  f(x-xQ,  y-yQ)  ^  F^fx»  fy)  exP  ^_i2lT^fxxo  +  fyyo^ 


The  triangle  function  tri(x),  tri(y)  is  defined  as: 


tri  (x) 


tri  (y) 


1  -  1*1 
10 

i  -  lyl 
lo 


W  <  i 

elsewhere 

|y|  U 

elsewhere 


Equation  (2)  for  the  Fourier  transform  and  its  associated  conditions 

for  the  pupil  function,  unit  amplitude  and  f-d  can  be  applied  to  the  two- 

dimensional  triangle  function  where  now  we  let  t(x  ,y  )  be  the  amplitude 

o  o 

transmittance. 


t(xo*yo)  “  trl  (xo>  tri  (yo> 

then  the  Fourier  transform  becomes 

A  exp  [i  —-(x^  +  y* )  ] 
g2(xf,yf) - Hf - 


F  {t<Vyo)} 


where: 


(6) 


(7) 


F  {t(x  ,y  )}  a  sinc^  (f  )  sinc^  (f  ) 


(8) 


and  again:  =  (j^) ;  f  =  (jj) 


are  the  spatial  frequencies,  f  =  focal  length  of  the  converging  lens; 

X  ■  wavelength  of  Helium-Neon  laser;  x^  and  y^  are  coordinate  positions 
in  the  Fourier  transform  plane.  Since  we  measure  (photograph)  the 
intensity  distribution  of  g2Cx£,y^)  we  can  neglect  the  quadratic  phase 
factor  that  appears  in  eq  (2)  and  (7),  The  resulting  intensity  for  the 
triangle  function  becomes 

XTri^Xf ,yf)  “  ^slnc2  <fx)  Sine2  (fy)}2  (9) 

C.  Lowpass  filter  "ringing" 


It  is  noticed  that  considerable  amount  of  "ringing"  and  some 
blurring  of  inverse  Fourier  transform  Images  occurs  when  sharp  lowpass 
frequency  cut-off  filters  are  used.  This  effect  is  seen  in  the  results 
of  the  Analogic  array  processor  IHFT  images.  We  too  were  able  to 
simulate  this  "ringing"  with  our  sharp  lowpass  optical  filters.  A 
sharp  lowpass  filter  (SLPF)  can  be  represented  by  a  transfer  function 
H(u,v)  in  the  Fourier  transform  plane,  where 


H(u,v) 


r(u,v)  <  rQ 
r(u,v)  >  rQ 


(10) 


rQ  is  the  radius  of  the  rectangle  function.  The  complete  function  is 
obtained  by  rotating  the  rect  function  about  its  axis.  This  generates  an 
"ideal"  lowpass  filter  with  sharp  edges.  The  radius  of  the  filter  is 


given  by 


r  -  [u2  +  v2]  **  (11) 

o 

This  circle  super-imposed  on  the  Fourier  transform  plane  (concentric 
with  it)  provides  the  interior  zone  for  passing  the  low  frequencies.  As 
rQ  is  increased  higher  frequencies  are  passed;  eventually  rQ  is  so  large 
that  all  frequencies  pass. 

The  "ringing"  phenomenon  can  be  discussed  in  terms  of  the  convolution 
theorem.  In  the  frequency  plane  we  have  that 

G(u,v)  »  H(u,v)  F  (u,v)  (12) 

where: 

H(u,v  )  »  transfer  function  and  F(u,v)  ■  Fourier  transform  of  the  Input 
function  f(x,y).  Using  the  convolution  theorem  one  gets 

g(x,y)  -  h(x,y)  *  f(x,y)  (13) 


where 


g(x,y)  ■  output  function  and  h(x,y)  ■  inverse  Fourier  transform  of 
the  transfer  function.  Since  the  transfer  function  is  the  rectangle 
function  its  inverse  is  the  sine  function.  Therefore,  the  convolution  of 
this  sine  function  with  f(x,y)  will  yield  a  series  of  rings  depending  on 


the  radius  of  h(u,v).  The  number  of  rings  produced  is  inversely  proportional 
to  the  radius  of  H(u,v).  For  lowpass  filtering  more  ringing  occurs  as  the 
radius  gets  smaller.  When  the  radius  of  H(u,v)  is  large  enough  to  include 
all  the  high  frequencies  the  ringing  disappears;  but  so  does  the  lowpass 
filtering  effect, 

D.  Elimination  of  "ringing" 

It  is  possible  to  lessen  or  remove  the  ringing  in  lowpass  filtering. 

The  method  is  basically  one  of  smoothing  the  edges  of  the  lowpass  filter. 
Special  filters  can  be  constructed  for  this  purpose.  We  shall  give  three 
examples  of  such  filters. 


where  the  order  of  the  transfer  function  H(u,v)  is  given  by  n;  rQ  is  the 
cut-off  frequency.  By  rotating  the  H(u,v)  function  about  its  axis  one  gets 
the  three-dimensional  form.. 
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III.  EXPERIMENTAL  RESULTS 

Shown  in  figures  (3  -  8)  are  the  results  of  the  comparison  between 
the  hardware  Analogic  array  processor  Fourier  transforms  and  those 
obtained  using  optical  Fourier  transform  methods.  In  order  to  better 
compare  the  two  methods  the  35  mm  Air  Force  transparancies  (negatives) 
were  converted  into  positives  as  ours  were  before  printing  the  hardcopy 
photographs.  In  this  way  the  backgrounds  for  both  cases  is  the  same  (i.e. 
dark) • 

Figure  (3)  shows  the  Air  Force  results  obtained  for  the  square  aperture 
input  using  the  Analogic  array  processor.  Figure  (3b)  is  the  HFT  (hard¬ 
ware  Fourier  transform)  of  the  square  aperture  in  (3a).  This  result  is 

consistent  with  its  optical  counterpart  OFT  (optical  Fourier  transform) 

» 

as  shown  in  fig.  (4b)  for  the  square  aperture  input  (4a).  Only  a  slight 
scale  size  difference  is  observed  due  to  the  smaller  square  used  in  (4a) . 

The  slightly  smaller  square  was  used  in  the  optical  set-up  in  order  to 
utilize  the  center  of  the  lens  and  thereby  eliminate  possible  aberrations. 
Figure  (3c)  shows  the  IHFT  (inverse  hardware  Fourier  transform)  of  the 
lowpass  filtered  Fourier  transform  shown  in  (3d).  The  optical  comparison 
is  shown  respectively  in  figure  (4c)  and  (4d).  Some  differences  are 
observed  in  these  two  methods  and  are  as  follows:  (1)  more  internal 
"ringing"  is  shown  in  IHFT  (3c)  than  the  IOFT  (4a).  This  could  be  due  to 
the  fact  that  when  photographing  (4c)  the  optical  IOFT  the  film  quickly 
saturates  in  intensity  at  the  center,  thereby  giving  the  "ringing"  less 
contrast.  However,  in  the  IHFT  the  displayed  results  were  presented  in 
the  logarithmic  scale  thereby  enhancing  the  "ringing".  Visually  we  observed 
the  same  "ringing"  phenomenon.  The  asymmetry  is  the  "ringing"  shown 
in  figure  (3c)  IHFT  was,  however,  not  observed.  We  believe  this  could  be 


due  either  to  a  hardware  or  software  problem.  Another  possibility  is 
when  the  square  aperture  was  digitized  for  input  to  the  Analogic  array 
processor  it  is  possible  the  illumination  of  the  35  mm  transparency  was  not 
uniform. 

The  IOFT  (inverse  optical  Fourier  transform)  of  the  highpass  filter 
shown  in  fig  (5d)  is  presented  in  (5c).  Again  "ringing"  and  harmonic 
mode  oscillations  of  the  "ringing"  are  observed  in  the  optical  case.  The 
hardware  IHFT  filter  shown  in  fig  (5a)  is  for  the  highpass  filter  in  (5b) . 
Although  some  edge  "ringing"  is  observed  we  do  not  see  the  same  Interior 
mode  "ringing"  observed  optically.  We  believe  that  either  fig  (5a)  is  not 
the  actual  IHFT  of  (5b)  or  the  exposure  of  (5a)  was  not  presented  on  a 
logarithmic  scale  as  was  done  for  fig  (3c)  in  the  lowpass  IHFT. 

The  six-pointed  star  aperture  and  its  hardware  Fourier  transform  are 
shown  in  fig  (6a)  and  (6b)  respectively.  The  lowpass  IHFT  of  the  star  is 
shown  in  fig  (6c) .  Again  one  observed  the  "ringing"  as  well  as  some 
asymmetry  in  the  "ringing"  in  the  inner-most  rings.  The  optical  results 
to  compare  against  the  star  are  shown  in  fig  (7) .  In  this  case  a  clear 
cut-out  of  a  hexagon  was  used  to  simulate  the  star  shwon  in  fig  (6a) .  This 
hexagon  fig  (7a)  via  the  Fourier  theorems  (discussed  earlier)  should 
replicate  to  a  good  approximation  fig  (6a) .  It  was  also  easier  to  fabricate 
accurately  than  a  six-pointed  star.  The  optical  Fourier  transform  of  (6a) 
is  shown  in  fig  (6b).  The  inverse  optical  Fourier  transform  for  the  low- 
pass  filter  is  given  in  fig  (7a).  "Ringing"  is  observed  as  before  and  due 
to  intensity  saturation  the  central  "ringing"  observed  visually  did  not 
photograph  well.  Again,  we  did  not  observe  any  asymmetry  in  the  central 
region  as  seen  on  the  logarithmic  scale  in  fig  (6c).  We  believe,  as  before, 
that  this  asymmetry  is  the  result  of  non-uniform  illumination;  hardware  or 
software  problems  in  the  array  processor  system* 


The  high  and  lowpass  opCical  filters  used  were  proportionaly  the 
same  size  as  those  used  in  the  hardware  analysis.  We  tried  to  use 
(roughly)  the  same  frequency  cut-offs  in  both  cases.  In  addition  we  tried 
to  eliminate  observed  "ringing"  by  using  smooth  high  and  lowpass  optical 
filters  instead  of  sharp  ones.  Figures  (8a)  and  (8c)  show  the  inverse 
Fourier  transforms  after  using  the  smooth  filters  on  (8b)  and  (8d) 
respectively.  One  observes  that  the  "ringing"  nearly  disappears.  These 
filters  were  produced  by  allowing  a  slight  blurring  of  the  boundary  between 
the  opaque  and  transparent  regions  on  the  glass  filter  slide.  The  results 
obtained  are  discussed  further  in  the  next  section. 
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IV.  DISCUSSION  AND  CONCLUSION 


The  comparison  of  the  optical  Fourier  transform  against  the  Analogic 
array  processor  Fourier  transform  of  the  square  and  star  aperture  showed 
that  both  similarities  and  some  differences  were  present.  We  shall  discuss 
each  separately  and  draw  some  conclusions. 

A.  Simularities 

The  Fourier  transforms  of  both  the  square  and  star  were  essentially 
the  same  for  the  HFT  and  OFT  methods.  There  are  advantages  in  the  hard¬ 
ware  method  when  displaying  the  logarithmic  of  the  intensity  rather  than 
a  photograph  of  the  Fourier/or  inverse/  transform.  Both  methods  also 
displayed  "ringing"  in  their  inverse  low  and  highpass  filtered  transforms. 
The  relative  size  of  the  low  and  highpass  filters  was  about  the  same  for 
both  methods. 

B.  Differences 


The  amount  of  "ringing"  varied  in  both  methods.  Here  the  "ringing" 
could  be  better  observed  when  displayed  by  the  hardware  (logarithmic) 
method.  Variations  in  the  amount  of  "ringing"  really  depend  on  the  exact 
frequency  cut-off  for  the  low  and  highpass  filters.  The  precise  cut-offs 
could  not  be  easily  matched  since  we  were  not  given  any  data  on  the 
hardware  parameters.  The  important  point  is  that  for  the  sharp  cut-off 
frequencies  "ringing"  was  observed  in  both  cases.  However,  an  important 
difference  is  the  observed  asymmetry  in  the  "ringing"  in  the  hardware 
method.  This  asymmetry  was  not  observed  at  all  in  the  optical  method. 


Some  conclusions  can  be  drawn  from  this  comparison.  The  "ringing" 
observed  is  to  be  expected  based  on  the  so-called  Gibbs  phenomenon.  The 
error  due  to  the  Gibbs  "ringing"  depends  on  where  one  truncates  the 
frequency  spectrum  before  performing  the  inverse  Fourier  transform.  The 
inclusion  of  more  high  order  frequencies  greatly  reduces  the  "ringing". 

When  performing  both  low  and  highpass  filtering  the  "ringing"  can  be 
minimized  by  using  a  "smooth"  filter  as  discussed  earlier  in  this  report. 
For  example,  an  array  processor  could  program  into  its  filtering  an 
exponential  decay  curve  which  would  pick  up  small  amounts  of  high  frequen¬ 
cies.  The  asymmetry  in  the  "ringing"  should  be  further  studied  to  deter¬ 
mine  its  exact  cause.  This  requires  the  array  processor  which  we  do  not 
have  in  our  laboratory.  There  are  pros  and  cons  for  each  method.  For 
example,  the  optical  Fourier  transform  due  to  its  parallel  processing  is 
faster  than  the  array  processor,  given  an  input  transparency  one  can 
readily  get  its  optical  F.  T.  using  a  rather  inexpensive  lens.  An  array 
processor  would  be  a  considerably  more  expensive  method  to  get  a  simple 
F.  T.  of  an  image.  On  the  other  hand,  manipulation  of  the  Fourier  trans¬ 
form  plane  with  various  Modulation  Transfer  Functions  (so  long  as  they 
can  be  obtained  analytically)  is  easier  (and  more  accurate)  to  perform 
with  the  array  processor.  Finally,  the  DIP  system  is  a  potentially  power¬ 
ful  and  important  digital  image  processing  system  capable  of  many  detailed 
Fourier  transform  studies.  It  should,  however,  be  further  compared  against 
optical  methods  as  indicated  by  this  preliminary  research  study. 
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Analogic  Array  Processor  Results 
(Air  Force  Hardware) 


Square  aperture  Input 
(Air  Force) 


(b) 

Fourier  transform  of 
square  aperture  (a) . 


Inverse  Fourier  transform 
of  lowpass  filtered  FT 
(b). 


Lowpass  filter  used 
on  Fourier  transform 
(b). 


Figure  (3) 


Optical  Processor  Results 
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Inverse  Fourier  transform  Lowpass  filter  used  on 
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Figure  (4) 


Analogic  Array  Processor  Results  (a.b) 


(Air  Force  Hardware) 
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Inverse  Fourier  transform 
of  highpass  filtered  FT 
Shown  in  Fig.  (3b). 


Highpass  filter  used 
on  Fourier  transform 
Fig.  (3b) 


Optical  Processor  Results  (c,d) 


(c)  (d) 

Inverse  Fourier  transform 

Highpass  filter  used 

of  highpass  filtered  FT 

on  Fourier  transform 

shown  in  Fig.  (4b). 

Fig.  (4b). 

Figure  (5) 
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Analogic  Array  Processor  Results 
(Air  Force  Hardware) 


(a) 

Star  aperture  Input 
(Air  Force) 


(b) 

Fourier  transform  of 
star  aperture  (a) 


Figure  (6) 


Optical  Processor  Results 


(a) 

Inverse  Fourier  transform 
of  smooth  lowpass 
filtered  FT  shown  in 
Fig.  (7b). 


Smooth  lowpass  filter 
used  on  FT  shown  in 
Fig  (7b)  for  simulated 
star  aperture. 


(c) 

Inverse  Fourier  transform 
of  smooth  highpass  filtered 
FT  shown  in  Fig.  (7b). 


(d) 

Smooth  highpass  filter 
used  on  FT  shown  in 
Fig  (7b)  for  simulated 
star  aperture. 


Figure  (8) 


